Influence of Landau-level mixing on Wigner crystallization in graphene 
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Graphene, with its massless hnearly-dispersing carriers, in the quantum HaU regime provides an 
instructive comparison with conventional two-dimensional (2D) systems in which carriers have a 
nonzero band mass and quadratic dispersion. We investigate the influence of Landau level mixing 
in graphene on Wigner crystal states in the n**^ Landau level obtained using single Landau level 
approximation. We show that the Landau level mixing does not qualitatively change the phase 
diagram as a function of partial filling factor v in the n^^ level. We find that the inter-Landau level 
mixing, quantified by relative occupations of the two Landau levels, pn+i/pn, oscillates around 2% 
and, in general, remains small (< 4%) irrespective of the Landau level index n. Our results show 
that the single Landau level approximation is applicable in high Landau levels, even though the 
energy gap between the adjacent Landau levels vanishes. 

PACS numbers: 73.20.Qt, 73.43.-f 

I. INTRODUCTION 

Wigner crystallization, where the density profile of carriers in a system develops a periodic spatial modulation 
spontaneously, is a classic example of interplay between (classical) repulsive potential energy and the (quantum) 
kinetic energy associated with localization of carriers as the density of carriers is varied. Although predicted in 
1934i this phenomenon has defied direct experimental observation in bulk systems and conventional 2D systems. 
In quantum Hall systems, where the kinetic energy of carriers is quantized and quenched, Wigner crystallization 
is induced by a competition between the electrostatic and exchange interactions as the partial filling factor v in 
a given Landau level is varied. (In the quantum Hall regime, Wigner crystallization depends only on the filling 
factor and can occur at any carrier densityi^) Wigner crystallization in the lowest Landau level has been inferred 
via transport measurements^ and the anisotropic transport observed^ in high Landau levels can be interpreted^ in 
terms of anisotropic Wigner crystal ground states. A direct observation of the Wigner crystal, via local carrier density, 
however, has not yet been possible. Graphene, with its massless carriers on the surface, is a unique and ideal candidate 
for this purposei^ Recent studies, using Hartree-Fock mean-field theory in the single-Landau-level approximation 
(SLLA)i2 or exact diagonalization in the single-Landau-level subspacc''^"'^ have predicted that Wigner crystal states 
will appear as ground states over a range of partial filling factor i/ in a given Landau level. In this paper, we examine 
the validity of the single-Landau-level approximation. 

Let us first recall the relevant results for a conventional 2D system in perpendicular magnetic field B with partial 
filling factor u < I'm the Landau level n. Thus, the actual filling factor for spinless carriers (with no other degeneracies) 
is n + v. For this system, the difference between energies of the adjacent Landau levels is AEn = -En+i — En = fiiOc 
where Uc — eB/mc is the cyclotron frequency, m ~ 0.5me — O.lme is the band mass of the carriers, and me is the 
bare electron mass. We remind the Reader that Ai?„ — /mlg is (approximately) the quantum kinetic energy of 
a particle with mass m in a box with size Lb = a/ hc/eB. The Coulomb interaction that causes transitions between 
different Landau levels has a typical energy scale Vc — e^/els where e ~ 10 is the dielectric constant. Therefore, 
the ratio of these two energy scales, Vc/AEn — Is/aB where qb — h'^e/me'^ is the Bohr radius of the carriers in the 
material. Since as is independent of the magnetic field and the magnetic length Is oc 1/^/B, as i? ^ oo the amplitude 
for inter-Landau level transitions vanishes and the SLLA becomes a good approximationii^ A corresponding analysis 
for graphene shows the stark difference bet ween the two systems. The gap between the adjacent Landau level energies 
in graphene is AEn = En+i — En = ftw[\/2(n -I- 1) — \/2ri\ where uj = vg/Ib is the cyclotron frequency, vg ^ c/300 
is the speed of massless carriers in graphene, and c is the speed of light. It follows that the ratio 

_ Vc ^ 1 
^" AEn ehvG ^2(n-f 1) - V2^ ' 

is independent of the magnetic field and diverges, gn ^ V2n, as n — s- oo. Therefore, inter-Landau level transitions 
become increasingly important as the Landau level index n increases, irrespective of the magnetic field; even in the 
lowest Landau level, the ratio g ~ e^/ehvG = etc ~ 1 is not small (ac is the fine structure constant for graphene). 
This analysis suggests that the SLLA is not reliable in graphene for any B and that it gets worse with increasing n 
since the energy gap AEn ~^ 0. In the following we show that, contrary to the expectations from a simple analysis 
presented above, the effect of Landau level mixing in graphene remains small and SLLA remains applicable. 
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The outline of the paper is as follows. In Sec. [Til we briefly describe the Hartree-Fock approximation with Landau- 
level mixing and outline our approach. The details presented in this section are essentially identical to those in our 
earlier work.— In Sec. Illli we present the results obtained without and with Landau-level mixing. We find that the 
Landau-level mixing does not qualitatively change the phase diagram of the system. We quantify the mixing using 
off-diagonal self-energy and relative occupation of Landau levels n and n + I. We compare the results for Landau 
level mixing as a function of n in graphcnc with those for conventional 2D systems. We summarize our conclusions 
in Sec. |TVl 



II. MICROSCOPIC HAMILTONIAN AND HARTREE-FOCK APPROXIMATION 



Let us consider graphene in a strong perpendicular magnetic field B in the quantum Hall regime. The single-particle 
states of the non-interacting system are given by \n, k, a) where (n, k) denote the Landau level and intra-Landau level 
indices, and cr = ± correspond to the two inequivalent valleys, K and K' = — K, in the Brillouin zone. The details 
presented in this section follow closely Ref.[l3]. The Hamiltonian for the system, including the Coulomb interaction 
is 



H = iv^^(£;„ - /i)p:^:^(o) + ^ E ^(q)-^«i.«2(-q)/5^1:nUq)-^"3.n4(-q)/5^r,«!(q) 



(2) 



where A is the area of the sample, /i is the chemical potential, V{q) ~ 27re^/eg is the Coulomb interaction in graphene 
(e ~ 2 - 5), and 



(3) 



k,k' 



with cjj^^(c„fe(j) representing the creation (annihilation) operator for state \n,k,a). Eq.Q is related to the density 
matrix operator in the momentum space 



(4) 



where the form factor for graphene (with n, n' > 0) is given bj*** 



(5) 



We recall that J-n^n' (q) is a linear combination of the form factors for a conventional 2D system 
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-Fn>«'(q) = V — r 



(iq-x - qy) 



V2 



(n—n') 



(6) 



where L™(x) is the generalized Laguerre polynomial and -F„<n'(q) = -F*, „(— q)- 

The derivation of the mean-field Hamiltonian using Hartree-Fock approximation is straightforwardi^ and gives 

an CTi 71 j q 

where C/crini,<T2n2 (q) = -ffcrini,<T2n2 (q) ^ -^crini ,cr2n2 (q) ■ Thc sclf-consistcnt elcctrostatic and exchange potentials are 
given by 



-ffCTini,o-2n2 (q) ~ ^o"i,(T2 ^ ^ ^nin3,n2n4{(l)Pn3,ni{ q) J 
n3n4cj 

^crini,cr2n2 (q) — ^ ] ^nin3.n2n4 (q) ['^o'l ,CT2 Pna .nli (~q) + ^cri ,CT2 Pns'.nli (~q)] ) 



(8) 
(9) 
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where 



-^nin3,n2n4 (q) 



(2^ 



(10) 
(11) 



and Pn\'^^2^(\) ~ {Pn\'^2^°L}) ^"^^ density matrix elements which should be determined self-consistently from Eq.(l7]). 
The density matrix is obtained from the equal-time limit (r — s- 0^) of the single-particle Green's function 



The equation of motion for the Green's function in Fourier space is given hy^ 

0-3 "3 q' 

and the Hartree-Fock self-energy matrix is (p = q — q') 

2^0-1 ni ,0-3713 (q, Q ) = ^ ^ { [-^nimi,n3m3 (~P)P'n3,mi (P) ^ ^"imi ,n3m3 (~P)Pm3 ,mi (P)] ^o'i,'T3 



(12) 
(13) 



— Xr, 



>-nimi,n3m3( P)Pm3 ,mi (P)'5o-3 ,01 } 6 ^ '^'^ i (14) 

where we have defined p„3^„n(p) = So-pf^g „j (p). 

In single Landau level approximation for the n**^ level, all Landau level indices in Eq. (ll3p are the same, ni — n2 — 
— n. To account for the inter-Landau level transitions, we restrict the indices to n and n+l. The Green's function 
in the Landau-level space then becomes a 2 x 2 matrix. 



'^n+l,n "-^n+l.n-l-l 



(q, «t^n) 



(15) 



and similarly the self-energy matrix Yicn,(T2 (Qi q') is a 2x2 matrix in the Landau-level space. The equation of motion 
for the Green's function, Eq. (fT5|) . becomes 



■q,o = N + ^A G"^'"Hq,i^„) - [^-1^-3 (q,q') + E6^,c,'Saua-] G'^-^'"Hq',iw„) 



(16) 



where the kinetic energy matrix in the Landau level space is diagonal, E — diag(i?„, i?„+i). We solve Eg.p^ by 
obtaining the eigenvalues and eigenvectors 



X! [^'^i,^3(q>q') + -E(5q,q'5<,i,03 v;3(q',fc) = wfcK-i(q, fc). 



(17) 



esq' 



Here ^^(q, k) = [V£j!'„(q, fc), Vjj*„^]^(q, A:)] is the eigenvector with eigenvalue uj^. We can construct the self-consistent 
mean-field Green's function using these eigenvectors^ 



G'^--^(q,zc.„) = E 



^.,(q,fc)V;UO,fc) 



which, in turn, leads to the self-consistent density matrix 



(18) 



(19) 



where f{x) — 0{~'x) denotes the Fermi function at zero temperature. The chemical potential /i is determined by the 
constraint that the total occupation in the two Landau levels is equal to the partial filling factor. 



(20) 



Using the self-consistent density matrix (|19p. we calculate the Hartree-Fock mean- field energy E^f for various trial 
lattice configurations to obtain the ground state crystal structure. 
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III. RESULTS 



We consider mean-field Wigner crystal lattices with two primitive lattice vectors ai = (a,6/2),a2 = (0,6) and 
define the lattice anisotropy as j — b/a. Note that the triangular lattice (7 = 2/V3 — 1.15) and quasi-striped states 
(7 0) are special cases of the general anisotropic lattice defined by these primitive vectors. The lattice constants a 
and b are determined by the constraint that a unit cell contains electrons, and are given by a = Ib^/'^.ttN^/v^ and 
b — aj. The reciprocal lattice vectors are Qmn = mhi + nh2 where bi = (27r/a)(l,0) and b2 — (27r/a)(— 1/2, 1/7) 
are the reciprocal lattice basis vectors. We determine the optimal lattice structure by choosing the 7 (0 < 7 < 2/^/3) 
and Ne that minimize the mean-field energy Ehf- In the following, we denote crystals with one electron per unit 
cell, TVe = 1, as Wigner crystals and those with Ne > 2 per unit cell as bubble crystals . ^'^'^'^ We first calculate the 
self-consistent density matrix without Landau level mixing, Pn^n^i^l) 7^ and Pn+Vni'i) = = Pn+Tn+ii^)- then 
use that matrix as the initial point for the density matrix with Landau-level mixing. 




FIG. 1: (Color Online) Ground state energy per particle, measured in units of e^/ds, for different crystal structures in the 
n = (left) and n — 3 (right) Landau levels in graphene. The top panel (a) shows results without Landau level mixing, whereas 
the bottom panel (b) shows results with mixing. We see that, in each case, the ground state energy is lowered due to Landau 
level mixing.— Note that the phase-diagram is qualitatively unchanged. 

Figured] shows mean- field energy per particle for various lattice structures as a function of 1/ for Landau level n = 
(left) and n = 3 (right). We note that the ground state energy for a given v is lowered by the Landau level mixing, 
as expected from perturbation theory^i^ For n = 0, we find that a triangular Wigner crystal is the mean-field ground 
state with or without Landau level mixing. For n = 3, we find that the triangular lattice remains a ground state for 
higher values of v when the inter-Landau level mixing is taken into account. Overall, the phase diagram of the system 
remains qualitatively unchanged. 

The most visible effect of inter-Landau level mixing is the systematic up-shift of critical values of v at which 
transitions from one crystal structure to another occur. For example, at n = 3 the transition from an isotropic 
Wigner crystal to an A^e = 2 anisotropic bubble state occurs at v ^ 0.20 without Landau-level mixing; this critical 
value is shifted upwards io v ^ 0.25 when the mixing is taken into account (Figure [T|). This shift is also visible in the 
lattice anisotropy 7(1^) for the ground state crystal structure, shown in Fig. [21 At small v, the lattice is triangular and 
7 = 2/\/3 — 1.15 is a constant. At higher values of the anisotropy increases leading to a quasi-striped structure for 
the ground state. We see from Fig. [3 that the region of stability of the triangular lattice increases when inter-Landau 
level transitions are taken into account. 

Results in Figs. [T] and [2] suggest that the effect of Landau-level mixing is not dominant in higher Landau levels, 
even though the energy gap between adjacent Landau levels becomes smaller. To understand this unexpected result, 
we recall that the inter-Landau level transitions from n — > n -|- 1 are determined by the off-diagonal self-energy matrix 
elements and the gap between adjacent Landau levels, Scm.o-n+i/AiJn. It follows from Eas. (|14ll0llip that for large n 

[n + l) Vn + 1 

because .7>i.„+i ~ l/VrT+T. We find that this asymptotic behavior is reproduced by our results. We quantify the 
Landau-level mixing by the ratio of relative occupations of the two levels in question, pn+i/ Pn where pm — X^o- Pm'ni(O)- 
Left panel in Fig. [3] shows the ratio Sn^n+i/ASn as a function of Landau level index n for graphene (solid red) and the 
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FIG. 2: (Color Online) Ground state lattice anisotropy 7(i^) in graphene for n = 2 (top) and n = 3 (bottom) Landau levels. 
The solid (red) line shows the result without mixing and the dashed (green) line indicates the result with Landau level mixing. 
7 = 2/%/3 = 1.15 corresponds to a triangular lattice whereas 7 < 0.5 corresponds to highly anisotropic Wigner crystal or 
quasi-striped states. 



conventional 2D system (dotted green) at partial filling factor v — 0.5. We see that the ratio Scrn,crn+i(0, Qoi)/Ai?„, 
for typical ofF-diagonal self-energy matrix element in graphene, is smaller than 4%. In contrast to this, the ratio 
and the self-energy for a conventional 2D system decreases monotonically, since Ai?„ = hujc is independent of n, 
and is well-described by a + 1 dependence at large n. We recall that this ratio for a conventional 2D system 

depends on the magnetic field B or the magnetic length Ib- Our results are for g — Ib/clb ~ 0.67 or Ib ^ 35 A. 
(This g = Vc/ for a conventional 2D system is equal to the g — ac in graphene with e — 3.3 as the dielectric 
constant.) The right panel in Fig. [3] shows the corresponding relative occupations for graphene (solid red) and the 
conventional 2D system (dotted green). The fact that this ratio, in the presence of inter-Landau level mixing, is small 
(pn+i/pn < 4%) provides complementary support for the validity of SLLA in graphene. 




FIG. 3: (Color Online) Left: E+,„;+,„+i/A£'„ as a function of n for partial filling factor v — 0.5 in graphene (solid red) and 
a conventional 2D system (dotted green). In graphene, the ratio is small for all n and shows that SLLA is applicable even 
in high Landau levels when AEn 0. In conventional 2D systems, the ratio decays monotonically as l/\/n + 1 for large n. 
Right: Corresponding ratio of occupation numbers pn+i/pn for graphene (solid red) and conventional 2D system (dotted green) 
provides further support for the validity of SLLA in high Landau levels. 



IV. DISCUSSION 



In this paper, we have investigated the effects of inter-Landau level transitions on Wigner crystal mean-field states 
in graphene obtained using single-Landau- level approximation.^" Our results show that the Landau- level mixing does 
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not qualitatively change the phase diagram of the system, although it shifts upwards the critical values of filling 
factor V at which transitions from one lattice structure to another occur. We quantify the Landau-level mixing in 
terms of off-diagonal self-energy and relative occupation numbers, and show that it remains small as a function of the 
Landau level index n. Thus we conclude that SLLA provides a reliable description of Wigner crystal ground states 
in graphene. 

We emphasize that our results for graphene are independent of the magnetic field B. For conventional 2D systems, 
the Landau-level mixing depends on the magnetic field and can be important at weak fields B < Be when the 
magnetic length becomes larger than the Bohr radius of the massive carriers, Ib ^ as for B < Be- The absence of 
a corresponding critical field Be in graphene is due to the massless nature of the carriers. Our conclusions do not 
depend, qualitatively, on the range of the interaction T^(q) because the large-g scattering is strongly suppressed by the 
form factors .F(q) that decay exponentially with q. In this paper, we have ignored transitions to next-higher Landau 
levels^ n n + k, because the amplitude for them vanishes rapidly: Y^n^n+k/ ^Enk ^ \/ [n + + k)\k'^ ^ as 

n — » oo for fc > 2. This estimate follows from an analysis similar to that for Eq. pTj) and the observation that, in 
graphene, AEnk — E„-^.k — En oc k/y/n for large n ^ k. Therefore, it is sufhcient to consider the Landau-level mixing 
only between adjacent levels. 

Since carriers in graphene are on the surface, in contrast to those in the conventional 2D system, it is an ideal 
candidate for direct observation of the local carrier density structure.^ Our results provide further support for the 
existence of triangular Wigner lattice as the ground state at small i> and anisotropic ground states in high Landau 
levels for v — > 1/2.^-'''^^ A direct observation of carrier density in graphene in the quantum Hall regime will verify (or 
falsify) our conclusions. 
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